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The possibility that the energy density of the Universe is dominated by a network of low-tension 
domain walls provides an alternative to the commonly discussed cosmological constant and scalar- 
field quintessence models of dark energy. We quantify the lower bound on the number density of 
the domain walls that follows from the observed near-isotropy of the cosmic microwave background 
radiation. This bound can be satisfied by a strongly frustrated domain wall network. No fine- 
tuning of the parameters of the underlying field theory model is required. We briefly outline the 
observational consequences of this model. 



I. INTRODUCTION 

In the last few years, there has been a growing amount 
of evidence that an unknown negative pressure compo- 
nent ("dark energy" ) accounts for 70 ± 10% of the energy 
density of the Universe. Arguably the most direct evi- 
dence comes from the luminosity-redshift relation of the 
Type la supernovae [1,2]. Additional evidence is pro- 
vided by a combination of data [3]. Both the analysis 
of the formation of the large scale structure [4] and the 
studies of the baryon fraction in galaxy clusters [5] , com- 
bined with the Big Bang Nucleosynthesis calculations [6] , 
suggest the matter content significantly below the crit- 
ical density. On the other hand, the cosmic microwave 
background radiation (CMBR) power spectrum suggests 
a flat Universe [7] . Therefore the difference must be made 
up by an additional non-clustering energy density com- 
ponent. The relationship between the pressure p and the 
energy density p of the dark energy component is usually 
parameterized by p — wp, where the equation of state 
w could be time dependent. Theoretical considerations 
prefer w > — 1. A good fit to the supernova data can be 
obtained assuming a constant or slowly changing equa- 
tion of state satisfying w < —0.5 [3] or w < —0.6 [8], 
depending on the details of the analysis. 

The precise nature of the dark energy component has 
been a subject of intense theoretical speculation. The 
most popular dark energy candidates include the cos- 
mological constant, or vacuum energy, and the so-called 
quintessence, a nearly spatially homogeneous but time- 
dependent scalar field (see, for example, Refs. [9-11].) 
Another very interesting possibility is the "solid" dark 
energy, originally suggested by Bucher and Spergel 
[12,13]. In this case, one postulates that the dark en- 
ergy component possesses resistance to pure shear de- 
formations, guaranteeing stability with respect to small 
perturbations. (The stability condition is nontrivial for 
substances with negative pressure. It is violated, for ex- 



ample, by a perfect fluid.) 

A possible microphysical origin for solid dark energy is 
a dense network of low-tension domain walls [13,14]. This 
is attractive for several reasons. First, domain walls are 
ubiquitous in field theory, inevitably appearing in models 
with spontaneously broken discrete symmetries [15]. Sec- 
ond, domain walls, and the solid dark energy in general, 
have been shown to be compatible with the observations 
of large scale structure [13,16]. Finally, a static wall net- 
work has an equation of state w = —2/3 [17], consistent 
with all observational data. Despite these appealing fea- 
tures, the idea that a domain wall network could play 
the role of dark energy has not received much attention 
in the literature. In this paper, we investigate certain 
aspects of this scenario. We show that the requirement 
of the isotropy of the CMBR yields a strong lower bound 
on the number density of the domain walls. We then con- 
sider the implications of this bound for the field-theoretic 
models of the domain walls. We argue that the domain 
wall network has to be strongly frustrated to satisfy the 
CMBR constraint. If this is the case, the equation of 
state of dark energy is expected to be close to —2/3. 

II. CMBR CONSTRAINTS 

It is well known that the Universe is spatially homo- 
geneous and isotropic on large scales. The best evidence 
for this is provided by the observed near-isotropy of the 
CMBR. If a network of domain walls is present, the 
CMBR photons would acquire different gravitational red- 
shifts (or blueshifts) depending on the direction of their 
arrival [18]. Let us estimate this effect. 

Until z ~ O(l), the contribution of the walls to the 
energy density was subdominant and hence they did not 
significantly affect the CMBR anisotropy. Most of the 
anisotropy due to domain walls would be built up dur- 
ing the relatively recent epoch, z < 1. This allows us to 
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simplify the problem. Consider an observer at rest in the 
comoving frame. (After subtracting the effects due to the 
peculiar motion of the Earth, the CMBR anisotropy mea- 
sured in actual experiments corresponds to what would 
be measured by such an observer.) Let R denote the co- 
moving distance corresponding to z = zq < 1, and con- 
sider the sphere S of radius R centered on the observer. 
We assume that the CMBR photons arriving at the sur- 
face of this sphere are isotropic. The anisotropy induced 
when the photons are travelling from the surface of S to 
the observer can be simply estimated, up to corrections 
suppressed by powers of zo, using Newtonian theory. In- 
deed, the results of the general theory of relativity on 
subhorizon scales with receding velocities v <C c must be 
equally well described using Newtonian gravity. Recall 
that in the Newtonian picture, the observer is taken to 
be at "the center of the Universe" . Homogeneous matter 
around him creates a gravitational field which grows with 
distance as 

S(r) = -G— (p + 3p)rr. (1) 

Here we are following the general principle that the 
source of Newtonian gravity is the combination (p + 3p) 
and not the density itself (see, for example, [19]). The 
redshifts of distant objects in this picture result from the 
Dopplcr shift due to the expansion of the Universe, and 
the effect of the gravitational field in Eq. (1). In Ap- 
pendix A, we discuss the calculation of these effects in 
a homogeneous Universe, and show that the Newtonian 
calculation is accurate up to and including the second or- 
der terms in the relative velocity of the emitter and the 
observer. 

If domain walls are present, an additional, anisotropic 
gravitational field will be induced on top of (1), leading 
to anisotropies in the CMBR temperature. Denoting the 
Newtonian potential of the wall network by $(r), we can 
write the observed temperature difference between points 
1 and 2 on the surface of the sphere S as 

^2l^P = *, ri )- 4 M, ,2, 

where r 12 are the vectors pointing from the center of the 
sphere to points 1 and 2 and f 12 = r 12 /R- 

The Newtonian potential at point x from a single pla- 
nar wall is given by [20] 

c/>(x; n,x ) = -2nG N a n • (x - x ), (3) 

where a is the wall tension, xo is a point on the wall 
and n is unit vector normal to the wall. The force on 
a matter particle from the wall is repulsive. The total 
potential due to an arbitrary network of planar walls is 
the sum $(x) = ^\ </>(x; n^Xj). Because the potential 
(3) grows linearly at large distances, the walls outside 
the sphere cannot be neglected in the calculation of $. 
On the other hand, it is clear that the contribution of 



the walls lying outside the event horizon has to vanish. 
To model this effect, we introduce a regulator which is 
spherically symmetric with respect to the observer and 
cuts off the effects of the walls beyond certain radius 
R > R. (Notice that such a regulator is also necessary 
in the usual Newtonian treatment of a matter dominated 
FRW cosmology.) The dependence on the regulator will 
disappear in the final results. 

To see how the CMBR anisotropy generated by the 
walls depends on their distribution in space, let us con- 
sider a "one-dimensional" toy model in which all the 
walls are planar and perpendicular to the x axis. The 
anisotropy, defined in this model as the temperature dif- 
ference for the photons arriving from the ±x directions 
and measured by an observer at x = 0, is 

A$ = $(i?) - $(-i?). (4) 

Let N denote the number of walls between x = and 
x = R. Then, we obtain 

= F R + 4ttNG n <t ^R- JjJ2 x ^J ' ( 5 ) 

where 

(oo OO \ 

- X>/(z*) + X>/(z-i) (6) 
i=l i=l J 

is the gravitational force exerted by the walls on a unit 
mass placed at the origin, and Xi denote the positions 
of the walls to the right (i > 0) or to the left (i < 0) 
of the observer. The regulator function / in Eq. (6) is 
even, f(x) = f(—x), and smoothly cuts off the force 
due to the walls beyond the observer's event horizon: 
f(x) = 1 for |a;| < R, f(x) -» for \x\ > R. Using Eq. (5) 
and the analogous expression for <!>(— i?), we can obtain 
simple analytic estimates of the observed anisotropy. In 
our estimates, we will always assume N » 1; clearly, this 
condition has to be satisfied for the Universe to be even 
approximately homogeneous on large scales. In fact, the 
CMBR bound derived below will require N to be large, 
so our analysis is self-consistent. 

If the walls form a perfectly regular lattice structure 
with period L « fl, the leading dependence on R in 
Eq. (4) cancels out and one finds [18] 

A$ - 2itG n <jL. (7) 

This is exactly what one expects on physical grounds. 
Indeed, on scales much larger than L, a regular wall 
structure behaves like uniform dark energy, with its New- 
tonian potential approaching a symmetric parabola; all 
deviations from the parabola occur because of the gran- 
ularity of the structure on scales ~ L. However, Eq. (7) 
only holds if the wall structure is perfectly regular on 



2 



all scales up to the present size of the horizon*. Both 
causality considerations and the inherent randomness of 
the system make it extremely hard to believe that such 
a regular structure can be realized physically. 

In realistic models, the main contribution to the 
anisotropy comes not from the granularity of the wall 
network, but from long-wavelength fluctuations of the 
effective average energy density of the walls due to devia- 
tions from perfect regularity. As an illustration, consider 
an example in which the walls are displaced from their 
lattice positions by random amounts \Sxi\ <~ L/2. In this 
case, we find 

A$ ~ 2irG N <r(LR) 1/2 , (8) 

a much stronger anisotropy than indicated by Eq. (7). 
Even this model, however, is hardly realistic, since it re- 
quires that the number of walls to the left and to the 
right of the observer be identical. Removing even a sin- 
gle wall on one side of the observer leads, on average, to 
an anisotropy 

A$ - 2itGn<jR ~ 2ttG N (tNL. (9) 

In the absence of long-distance correlations between wall 
positions, the average difference in the number of walls 
to the left and to the right of the observer is of order 
\/N, leading to 

A$ ~ 2nG N aR 3 / 2 L- 1 / 2 ~ 2nG N aN 3 / 2 L. (10) 

The preceding discussion can be summarized by writ- 
ing the anisotropy in the form 

A$ ~ 2ttG N (tN v L, (11) 

where the value of the exponent v depends on the details 
of the wall configuration. What values of v correspond to 
realistic domain wall networks? To answer this question, 
let us consider a field theory with N v distinct vacua. Do- 
main walls form during a cosmological phase transition 
due to the fact that the field may choose different vacua 
in causally disconnected regions of space. Immediately 
after the phase transition, each causally connected re- 
gion contains either no walls, with probability p ~ 1/N V , 
or a single wall with probability 1 — p. First, consider 
the case when the wall network does not evolve (apart 
from conformal Hubble stretching) after the transition. 
In this case, the CMBR anisotropy induced by the walls 
today can be estimated as 

A$ ~ 2irG N <TN i ' 2 L^p{l-p). (12) 



'Although the final result is regulator-independent, the can- 
cellation of the independent terms requires a conspiracy be- 
tween the "nearby" walls at x < R and the "distant" walls at 
x ~ R. 



Thus, for moderate p, a non-evolving domain wall net- 
work induces the anisotropy of the size indicated by 
Eq. (10). The same estimate applies if the network does 
evolve, but the evolution does not make the structure 
more regular. On the other hand, if more regular wall 
configurations are favored dynamically and emerge as 
the network is evolving, the induced anisotropy could be 
weaker. While we are not aware of any numerical simu- 
lations that conclusively demonstrate such behavior in a 
given model, this remains a logical possibility. A reason- 
able lower bound on the induced anisotropy is provided 
by Eq. (9), since even a single "defect" in the regular 
wall structure inside the present Hubble volume would 
induce an anisotropy of that size. Below, we will derive 
the bounds on the parameters of the model using Eq. (12) 
(the "non-evolving network" case, v = 3/2) and Eq. (9) 
(the "regular network" case, v = 1). The constraints in 
any realistic model are expected to lie between these two 
limiting cases. 

While we have used a one-dimensional toy model to de- 
rive Eqs. (7 — 12), it is possible to show that the same es- 
timates hold for three-dimensional domain wall networks 
with the corresponding regularity properties. Thus, the 
temperature anisotropy created by domain walls inside a 
sphere of radius R for the non-evolving network case can 
be written as 

^ = 2™ G N <JN*' 2 L y/rfX=pj, (13) 

where a is a numerical coefficient of order unity that de- 
pends on the detailed properties of the network, L is the 
average separation between the walls, and N — R/L is 
the average number of walls crossed by a CMBR photon 
traveling from the surface of the sphere to the observer. 
For a regular network, we obtain 

^ = 2-kclG n cjNL. (14) 

This anisotropy is generated by the fluctuations of 
the wall number density at large (Hubble) distance 
scales. Predicting the power spectrum of the wall- 
induced anisotropy, as well as possible deviations from 
gaussianity, would require detailed knowledge of the net- 
work geometry. 

To be a viable dark energy candidate, the network 
of domain walls must have the average energy density 
p = 3cr /L = Ct w pcrit, where p cr it = 3Hq/8ttGn is the 
current critical density, and 0.6 < Q w < 0.8. (We as- 
sume that the walls are the only form of dark energy 
present: for example, the cosmological constant van- 
ishes.) To avoid conflict with the precise measurements of 
the CMBR anisotropy [21], the wall contribution should 
be at most at the level ST/T ~ 10~ 6 . Using Eq. (11), we 
find 

N> (2 x 10 5 a& 2 ) 1/(2_I/) 5 (15) 
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where b is the radius of the sphere in Hubble units, b — 
RHq w zq. In the case of a non-evolving network, we 
obtain 



iV>3x 10 10 a 2 6 4 p(l -p), 



< 



2Xl0 " 6 Ge V 3, 



a 2 b 3 p(l — p) 
0.15 



pc, 



a 2 b 3 p(l — p) 
while in the case of a regular network 
N>2x 10 5 ab 2 , 



(? ^ — 4 x 10 
ab 



-li 



GeV 3 , 



L<-3xl0 4 pc. 
ab 



(16) 
(17) 

(18) 

(19) 
(20) 

(21) 



Since the wall energy density becomes comparable to that 
of matter at z ~ 0.5, we will use the value b — 0.5 for our 
numerical estimates. The other two parameters entering 
the above bounds, a and p, depend on many factors, 
such as the geometry of the wall network, the properties 
of the underlying field theory model, etc. We will keep 
the dependence on these parameters explicit throughout 
the discussion. 

Our calculation has only included the anisotropy due 
to the walls at low redshifts, where Newtonian approxi- 
mation is applicable. The effects of walls at higher red- 
shifts could be included using the full general relativis- 
tic solution for the metric perturbation created by the 
walls in the FRW Universe [22]. In this formalism, the 
rcdshift is computed using the well known Sachs- Wolfe 
formula [23]: 



f) 

1 / sw 



$1° 



(22) 



where h^ v is the metric perturbation due to the walls 
< f ) = ftoo/2 is the "conformal" Newtonian potential, and 
x^ = (const+£, £e) is the unperturbed photon path. 
The subscripts e and o refer to the "emitter" (i.e., a 
point on the surface of last scattering) and the observer, 
respectively. Note that in an expanding Universe, even 
static domain walls induce a time-dependent metric per- 
turbation, and the third term in Eq. (22) does not van- 
ish. It was shown in [22] that in matter-dominated 



''"Explicitly, ds 2 = a 2 {rj) (gffi + h^dx^dx" , where rj is con- 
formal time. 



Universe, the anisotropy due to a wall at a redshift z w 
scales as (1 + z w )~ 3 / 2 (the redshift of a wall is defined 
as the redshift corresponding to the point on the wall 
closest to the observer.) Moreover, a constant average 
number density of walls in comoving coordinates im- 
plies that dN/dz w oc (1 + z w )~ 3 / 2 . Both these effects 
lead to a severe suppression of the effects of the walls 
at high redshifts. We have checked that the calculation 
of Ref. [22] leads to the same order-of-magnitude esti- 
mates (9) and (13) for the anisotropy created by the walls 
as the Newtonian treatment of our paper, justifying the 
latter. 



III. IMPLICATIONS 

Let us discuss the implications of the CMBR 
bounds (16 — 21) for the field-theoretic models respon- 
sible for the walls. Domain walls necessarily appear in 
theories with spontaneously broken discrete symmetries. 
If all the dimensionless parameters of the model are of or- 
der one, the wall tension is determined by the symmetry 
breaking scale v. The bounds (17) and (20) then imply 



< 



10 



a 2 /3pV3(l_p)l/3 



kcV 



and 



< 



4 x 10 2 



kcV 



(23) 



(24) 



for the non-evolving and the regular network cases, re- 
spectively. Note that the dependence on a and p is rather 
mild. In supersymmetric models, such low energy scales 
can be generated naturally and be radiatively stable, pro- 
vided that the SUSY breaking is communicated to the 
field(s) responsible for domain walls only by gravitational 
interactions. In this case, the natural value of the dis- 
crete symmetry breaking scale is given by v ~ F/Mpi, 
where VF is the scale at which SUSY is broken. The 
constraints (23) and (24) are satisfied if y/F < 10 4 - 10 5 
TeV, which is allowed phcnomcnologically if the breaking 
is mediated to the visible sector (Standard Model fields 
and their superpartners) by gauge interactions [24]. In 
this respect, the domain wall models of dark energy are 
much more attractive than the scalar-field quintessence 
models which contain a superlight scalar field with a mass 
of order 10~ 33 cV. In the latter case, it is difficult to un- 
derstand how such a low energy scale can arise from par- 
ticle physics and not be destabilized by radiative correc- 
tions. Moreover, the superlight scalar will in general me- 
diate a phenomenologically problematic new long-range 
force [25]. 

The simplest field theory model in which domain walls 
arise contains a single real scalar field with a Z„-invariant 
potential. (A well-known example is an axion of Peccei- 
Quinn models [26].) Numerical simulations [27,28] and 
analytic calculations [29] show that the wall network of 
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this model enters the so-called scaling regime shortly after 
its formation. In this regime, there is only one wall (or at 
most a few walls) per horizon volume at any given time. 
This is in gross contradiction with the bound (16) and, 
more generally, with the whole concept of the isotropic 
FRW cosmology. The problem could be avoided by intro- 
ducing a very large number (~ 10 8 ) of independent scalar 
fields, or considering models with a very large number of 
vacua. (In the context of the axion, the latter possibil- 
ity would correspond to a humongous value of the color 
anomaly of the Peccei-Quinn symmetry, N P q <~ 10 s .) 
Needless to say, both possibilities are extremely unattrac- 
tive. 

The scaling evolution of a domain wall network relics 
on the fact that the network can disentangle as fast as 
allowed by causality. In models with more complicated 
vacuum structure, the disentanglement process could be 
slowed down, since domains of the same vacuum would 
generally be well separated from each other in space. 
While no convincing simulations of domain wall networks 
exhibiting such behavior exist at present, a similar phe- 
nomenon has indeed been observed in the simulations 
of cosmic string networks [30]. Alternatively, the walls 
could experience strong friction forces, for example due 
to their interaction with the particles of dark matter [31]. 
In both cases, one expects substantial deviations from the 
scaling law. The networks of this kind are referred to as 
frustrated. For frustrated networks, the number of walls 
in the present horizon volume depends on their evolu- 
tion, and predicting it would require detailed numerical 
simulations. On the other hand, one can obtain a simple 
analytical upper bound on this number as a function of 
the wall formation temperature Tj. Assuming that the 
walls form during the radiation dominated epoch, and 
that the fields of the hidden sector responsible for the 
walls are in thermal equilibrium at approximately the 
same temperature as the visible sector fields' 1 ', we obtain 



T 2 



M Pl H 



(25) 



where z(T) is the rcdshift corresponding to the given 
temperature. (This bound is saturated by a non-evolving 
network, which in the present context can also be termed 
maximally frustrated.) In a model with no unnatural 
dimcnsionlcss parameters, Tf ~ v, and Eq. (23) leads 
to N h < 10 5 a~ 2 / 3 p~ 1 / 3 (l -p) -1 / 3 for the non-evolving 
case. This value does not contradict the CMBR con- 
straint (16) for reasonable values of a and p (for example, 
a ~ p ~ 0.1.) For the case of a regular network, Eqs. (24) 



*An upper bound on the temperature of the hidden sector 
fields is provided by Big Bang Nucleosynthesis. Unless the 
hidden sector possesses a large number of degrees of freedom, 
this bound can be satisfied with T bi d somewhat lower than, 
but of the same order as, T v i s . 



and (25) imply that N h <6- 10 6 a" 1 / 3 , although the true 
bound is probably somewhat lower since some of the walls 
are likely to be destroyed during the evolution leading to 
a regular structure. In any case, the bound seems to be 
compatible with the CMBR constraint (19). Thus, we 
conclude that in the presence of frustration, it should 
be possible to build realistic models of domain wall dark 
energy without any fine-tuning (apart from the tuning 
required to cancel the cosmological constant.) It would 
be very interesting to find explicit field theory models 
leading to a frustrated wall network. 

The upper bounds on N derived in the previous para- 
graph and the lower bounds in Eqs. (16) and (19) have 
to be nearly saturated to be compatible with each other 
for reasonable values of a and p. This observation leads 
to two interesting predictions. First, the wall network 
has to be close to maximal frustration, or, equivalcntly, 
be nearly static. Second, the CMBR anisotropy induced 
by the walls should be close to the current bounds, and 
therefore improved measurements of the anisotropy have 
a good chance of detecting the wall contribution if this 
model is realized. 



IV. EQUATION OF STATE 

It is well known that the equation of state of a maxi- 
mally frustrated (static) network of planar domain walls 
is w — —2/3. Above, we have argued that a domain 
wall network has to be close to the maximally frustrated 
regime to play the role of dark energy and be consistent 
with the CMBR anisotropy limits. One may worry, how- 
ever, that even small deviations from this regime may 
cause a large change in the value of w, making it an un- 
acceptable dark energy candidate. This turns out not to 
be the case. Let us demonstrate that the wall equation 
of state lies in the allowed range —2/3 < w < —1/2, 
regardless of the evolution. 

Energy conservation for domain walls reads d{pV) = 
—pdV — SE, where 5E is the energy lost by the network 
in the process of evolution. (The energy can be radiated 
away in form of gravitons, elementary excitations of the 
fields responsible for the walls, etc.) Let us tentatively set 
SE = 0. The wall equation of state can then be found 
from the dependence of the energy density of the wall 
network p w on the scale factor of the Universe a, 



w 



_ 1 d In p u 



3 din a 



(26) 



The number of walls in a fixed comoving volume scales 
as rj~ a , where rj is conformal time, dn = dt/a(t). When 
the wall network is static, a = 0. On the other hand, 
by causality there must be at least one wall per horizon 
volume at any given time, implying a < 1. The phys- 
ical density of the wall rest energy therefore scales as 
oc a~ 1 ri~ a . Neglecting the kinetic energy of the walls, we 
obtain 
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w 



2 a dlmj 

3 3 din a 



(27) 



During the matter dominated epoch the total pressure of 
the Universe p to t — 0, while after the wall energy takes 
over p^t becomes negative. It can then be shown that 
in an expanding Universe < din r]/ din a < 1/2, and 
therefore 



2 1 
- <w< --. 

3 ~ ~ 2 



(28) 



The lower limit is achieved when the domain wall network 
is static, while the upper limit corresponds to the scaling 
regime. Notice that even in the scaling regime, while the 
Universe is matter dominated, the wall energy density 
varies as oc a~ 3//2 . Hence, regardless of the details of the 
evolution, the wall network eventually comes to dominate 
the Universe and its equation of state always lies within 
the range allowed by the analysis of [3]. Of course, as al- 
ready mentioned, the network must be nearly maximally 
frustrated to satisfy the isotropy constraint, so we gener- 
ally expect w ~ — 2/3. In this case, the more restrictive 
bound obtained in [8] is also satisfied 5 . 

In the above derivation, we have neglected the en- 
ergy loss by the wall network SE and the kinetic energy 
of the walls Skin- Obviously, both approximations are 
valid for static networks. Moreover, numerical simula- 
tions [27,28] show that even in the scaling regime SE ~ 
and -Ekin <C -Brest throughout the evolution. This jus- 
tifies our assumptions. (Notice also that the condition 
SE = can be relaxed; the effect of finite SE/Sa is to 
shift w closer to -2/3, and Eq. (28) still holds.) On the 
other hand, the condition p to t < is essential for deriv- 
ing Eq. (28). For example, in the radiation dominated 
Universe the wall equation of state is w = — 1/3 in the 
scaling regime. 



V. CONCLUSIONS 

We have considered the constraints on domain wall 
models of dark energy from the observed near-isotropy 
of the CMBR. We have shown that these constraints can 
be satisfied by a strongly frustrated domain wall net- 
work. The scale of spontaneous symmetry breaking re- 
sponsible for the walls is expected to lie in the 10-100 keV 
range, and can arise naturally in supersymmetric theo- 
ries. This makes these models quite attractive from the 
particle physics point of view. 

Domain wall models of dark energy have important 
observational predictions. The dark energy equation of 
state is predicted to be close to —2/3. This value can 



5 Recent studies disfavor w = 
level [34] 



—2/3 at 95% confidence 



be clearly distinguished from the case of the cosmologi- 
cal constant by the SNAP experiment [32]. The CMBR 
anisotropy induced by the walls is likely to be close to 
the current bounds, and could be observable in the near 
future. In general, inhomogeneities of the wall distribu- 
tion are also expected to induce peculiar velocity flows 
on large scales. However, this effect is small: a network 
satisfying the CMBR constraints will produce velocities 
of order 300 m/s, which is not in conflict with current 
observations. 

The most important outstanding issue in making do- 
main wall dark energy models fully realistic is finding ex- 
plicit field theories which lead to highly frustrated wall 
networks. Frustration could arise as a result of the com- 
plex dynamics of the system. This possibility can only 
be addressed by numerical simulations. Such simulations 
are also necessary to make more detailed predictions of 
the phenomenological signatures of the walls, such as the 
power spectrum of the wall-induced CMBR anisotropies. 
At this time, we are aware of only one numerical study of 
frustrated networks [33] , whose usefulness is severely lim- 
ited by its insufficient dynamical range. Clearly, further 
work in this direction is necessary. An interesting al- 
ternative possibility is to consider domain wall networks 
whose evolution is slowed down by their interaction with 
dark matter. This idea was introduced in [31]. In the 
specific model of Ref. [31], the average wall velocity is 
determined by the ratio of wall and dark matter energy 
densities. In our case, this ratio is of order one, and the 
mechanism does not work. However, this class of ideas 
certainly deserves further investigation. 
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APPENDIX A: NEWTONIAN TREATMENT OF 
THE PHOTON REDSHIFT 

The rigorous way of computing the temperature 
anisotropy acquired by the photons on the way from the 
surface of last scattering to the observer is by using the 
techiques of general relativity On the other hand, the 
anisotropy induced when the photons are travelling in- 
side a sphere of a radius corresponding to z < 1 and 
centered on the observer can be simply estimated using 
Newtonian theory. (The estimate is expected to be ac- 
curate up to corrections suppressed by powers of z ) To 
illustrate how the Newtonian theory can be applied to 



G 



give the correct redshifts, accurate to the second order 
in v/c, let us explicitly compare the redshifts of a pho- 
ton emitted at some distance R computed by the two 
methods in the simplest case of exactly homogeneous, 
uniformly expanding Universe. 

From the point of view of general relativity, the ratio 
of the emitted and observed frequencies, ui/u>2 is simply 
the relative change of the scale factor of the Universe a 
during the time between emission and absorption, 



(Al) 



If the Universe is filled with a component at critial 
density with the equation of state w, the scale factor 
has a power law dependence on time, a <~ t n , where 
n = 2/[3(l + w)]. The velocity with which the emit- 
ter is receding from the observer at the time of emission 
ti is 



a(ti) h 



(A2) 



The time it takes for the light to travel to the observer, 
At = t 2 — ti, can be found from the equation 



rt-2 
Jti 



dt 



R 



l tl a{t) a(h)' 
Upon integration, we find 



R — n t™ — t\) 



1-n' 



so that 



t 2 



(1 - n)R + ti 



l/(l-n) 



(A3) 



(A4) 



(A5) 



Using this result, we find for the ratio of the scale factors 



(l-n)| + l 



n/(l-n) 



or, in terms of the recess velocity v — nR/t\, 



1 + v 



It \ \ n/(l—n) 

(1 — n) \ 



(A6) 



(A7) 



Finally, the ratio of the emitted to absorbed frequency is 



w(ti) 



1-v 



(n-1) 



n/{n-\) 



~ 1 - v + —v 2 + 



2n 



3(1 + w) 2 
1 - v+ -— -v 2 + ... 



(A8) 



Now let us compute the frequency redshift using New- 
tonian gravity. For that we take the observer to be at 
"the center of the Universe" , with the surrounding mat- 
ter inducing a gravitational field 



4tt 

g(r) = -Gy(p + 3p)rf. 



(A9) 



The frequency shift has two components, the kinematic 
Doppler shift and the blueshift because the photon falls 
into the potential well. The Doppler shift is simply given 

by 



LJ 2 



Doppler 



^i(l - v) 



~lu 1 (1-v + v 2 /2 + ...) (A10) 



The gravitational blueshift is computed as follows 

(w 2 -^i) Grav 



= A(f> = — I g{r)dr 
Jo 

G Y ( P + 3 P )-. 



(All) 



Since in the flat Universe the recession velocity obeys 
v(R) 2 /2 = GM(R)/R = 4ttG P R 2 /3, Eq. (All) could be 
rewritten as 



(, i , i ^Gra 
(L02 — Wl ) 



(1 + 3 W ) T - 



(A12) 



We observe that, by combining Eqs. (A10) and (A12), 
we recover precisely the result of Eq. (A8). Thus, in 
a homogeneous Universe, the Newtonian calculation of 
the photon redshifts due to the expansion of the Uni- 
verse is correct up to and including the second order 
terms in v. In Section II, we use the Newtonian pic- 
ture to estimate the anisotropics due to the presense of 
domain walls. Since the walls only become dominant for 
z ~ 1, this approximation should be sufficiently good for 
an order-of-magnitude estimate. This intuition is con- 
firmed by studying the full general relativistic calculation 
ofRef. [22]. 
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